%To interoplate and fft the tsunami v elocities
%lates date 19 June 2006

load c:\manoj\projects\tsunami\NewVel;
[OX,OY] = meshgrid([30.5:1:119.5]',[-29.5:1:28.5]'); %why these values ? = see the limit of lat1 and lon1
%Problem while interpolating: bilinear interpolation reduces the amplitude
%(probably only the shorter scale stuffs)

for i = 1:600,
    [IU_n(:,:,i)] = interp2(lon1',lat1',squeeze(IU(:,:,i))',OX,OY,'nearest');
end;
clear IU;
for i = 1:600,
    [IV_n(:,:,i)] = interp2(lon1',lat1',squeeze(IV(:,:,i))',OX,OY,'nearest');
end;
clear IV;